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Abstract 

We present results for higher order perturbative corrections to Compton scattering in the gener- 
alized Bjorken kinematics. The approach we have used is based on the combination of two tech- 
niques: conformal operator product expansion on the one side, and resummation of the fermion 
vacuum insertions with consequent restoration of the full QCD /3-function via the naive non- 
abelianization assumption, on the other. These are terms which are lost in the former approach. 
Due to the presence of the infrared renormalon poles in the Borel transform of the resummed am- 
plitude the latter suffers from ambiguities which reflect the asymptotic character of perturbation 
series. The residues of these IR renormalon poles give an estimate for the size of power corrections 
in deeply virtual Compton scattering. 
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1 Introduction. 



A measurement of deeply virtual Compton scattering (DVCS) [0], ||, j|, U should yield important 
information on the leading-twist non-forward parton densities^ which contain new information on 
the strong interaction dynamics and could open a new window for the exploration of the internal 
structure of the nucleon. To attain this goal one should have theoretical control of higher orders of 
perturbation theory as well as higher-twist corrections to the amplitudes. While QCD perturbation 
theory is well established up to the multiloop (2,3,4-loop) level, in practice all-order results for 
particular quantities are inaccessible due to the essential complexity of many loop calculations. 
The theory of power corrections is not settled completely. Even if Operator Product Expansion 
(OPE) or its generalizations — factorization theorems^ — exist for a particular process such that 
the power suppressed contributions can be expressed in terms of some multiparton correlators, 
the latter cannot be evaluated quantitatively mainly due to lack of a feasible non-perturbative 
approach. Therefore, in fact already a rough order of magnitude estimate for the power corrections 
would be extremely valuable. Furthermore the planned experiments on DVCS will obtain primarily 
data at rather limited Q 2 (this is especially true for CEBAF experiments). Thus there is legitimate 
concern that large higher-twist corrections might completely obscure the interpretation of such 
experiments at least for some kinematic regions. 

Recently, it has been realized that both of the above issues can be addressed by the study of 
the perturbative corrections in the large- Nf limit of QCD || [?J to the leading twist contributions. 
First of all, these calculations can be performed exactly due to the relatively easy algebra. Sec- 
ond, when supplemented by some assumptions, to be discussed below, it gives reasonable results 
when compared to exact perturbative quantities available. Last but not the least, by studying 
the ambiguities of the QCD perturbation series one can get some insight into the size of power 
suppressed corrections by the simple reason that only the sum of large orders of perturbation 
theory and higher twist contributions is free from ambiguities and thus physically relevant. 

One remark should be added to all that has been said above. Since the resummation of the 
vacuum insertions roughly corresponds to taking into account only effects related to the running 
of the coupling, all other radiative corrections are discarded at the same time. This turns out 
sometimes to be an unreliable approximation. On the other hand, by limiting oneself to the case 
when the /3-functions is zero, the so called hypothetical conformal limit, one has the advantage of 
conformal covarianceQ of the theory and can make use of conformal operator product expansion 

2 The naming "densities" which is used by us sometimes in the paper is not completely correct since it implies a 
probabilistic interpretation for the corresponding entries, however, the latter is lost for the non-forward functions. 
Therefore, it should be understood only as a synonym. 

3 The proofs of factorization for DVCS and diffractive meson production were given in Refs. ^). 

4 Of course, even for (3 = conformal invariance is broken by the renormalization of the field operators, however, 
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that allows to get strong restrictions on the form of the off-forward part of the massless amplitudes 
in higher orders of QCD perturbation theory. Therefore, below, when studying the higher order 
perturbative corrections we will combine both approaches, i.e., renormalon chain resummation 
and conformal OPE, which will amplify each other. 

The main objects of the present investigation are structure functions similar to those measured 
in deep inelastic scattering (DIS) but generalized to non-forward kinematics. They appear as the 
coefficients in the decomposition of the correlation function of two electromagnetic currents into 
independent Lorentz tensor structures which individually respect gauge invariance and are free of 
fictitious kinematical singularities 



T^(cu,C,Q 2 ) = i J d A ze iqz {h'\T {J^)J u (z)}\h) 



+ ^e^ p F>,C,Q 2 ). (1) 

The form factors F f (uj,() introduced above differ only by an overall constant from the ones 
considered by us in Ref. |J, namely F r = — \T r . The latter can be written in leading twist 
approximation as a convolution of the non-perturbative non-forward distributions and the pertur- 
batively calculable coefficient functions 

T f (ujX) = fdx £ i r r (w ) x,C,0 a |a.) i O r (a: > 0- (2) 

i=Q,G 

The former are defined as light-cone Fourier transformations of non-local string operators sand- 
wiched between appropriate hadronic states l O r (X,fi) ||: 

i O r (X,fi) = (/i'|0*(/in)r$[//n,An]0i(An)|/i) (3) 
= J dxe^ {x -°- iXx { i O r (x, ()6{x)6(l - x) - I O r (C - x, ()6(( - x)6{x + 1 - C)} • 

The limits of integration on the right hand side of Eq. (|j) can be deduced by studying the support 
properties of the function introduced above with Jaffe's approach ]TD| (see also [IT] for a recent 



discussion) which results^ in — 1 + ( < x < 1 [^]. Here i = Q,G runs over the parton species, 

r corresponds to different Dirac or Lorentz structures, depending on the spin of the constituents 

involved and $ is a path ordered exponential. The parton momentum fractions x and x — ( 

are the Fourier conjugated variables of the light-cone positions A and \l. The parameter ( is 

one can redefine the scale dimensions of the fields g] and embed them into the original conformal representation 
so that conformal covariance is preserved. 

5 In the original derivation the author of Ref. p] has used the perturbative analyticity approach for studying 

support properties of (multi)parton distributions jl^] (see also fl3f). 
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the skewedness of the distribution defined as a +-component of the t-channel momentum. The 
perturbative expansion of the coefficient function is given by 

^(u, x, C, Q 2 \a s ) = *T (0) (a;, x) + l T[ l} (u, x, (, Q 2 \a s ) + 0(a 2 s ), 

with the leading order (LO) hand-bag contribution^ Q T^(u,x) = uj/(xuj — 1) ± (x — > ( — x) 
with "+"-sign corresponding to unpolarized scattering and "— " to the spin- dependent case (and 
g T(q)(uj,x) = 0). We have followed above the conventions introduced by Radyushkin for DVCS 
[[J, namely, to = —2(pQ)/Q 2 , and A + = q + — Q + = (, with Q (p) and q (p') being the incoming 
and outgoing momenta of the photon (proton). They are related to those used by us in Ref. 

(which are closely connected to the variables adopted by the authors 01]) according to 
co = ~(PQ)/Q 2 = oo(2 - C)/(2 - O, V = C/(2 -Q,t= (2x - 0/(2 - 0, where the averaged 
momenta are introduced as follows P = p + p' , Q = \{Q + q) and A + = r]P + . An advantage 
of the first conventions is that the variable x acquires a simple partonic interpretation as the 
momentum fraction of the incoming constituent, while the variable t cannot be interpreted in this 
way. However, its range does not depend on the longitudinal asymmetry parameter rj, contrary 
to x. 

Considerable theoretical efforts have been undertaken recently to explore the properties of 
DVCS in leading order: The evolution kernels which govern logarithmic scaling violation have 
been evaluated in Refs. JTS], [H], [T], 0, ^, [H], [TH], [T^, [|. The solution of the renormalization 
group equation for the non- forward distributions given above was found in 0, [L4|], while explicit 
numerical studies were performed in |T7| , |14] , |20|| . However, the values of Q 2 for which the hand 
bag approximation can be trusted and gives satisfactory results need not be the same as in for 



usual forward DIS [21]. Therefore, this forces us to study the effect of higher twists and large 
order perturbative corrections to the amplitudes. 

The first step for any analysis beyond LO is the evaluation of the NLO coefficient functions. 
This issue was addressed first in [^2| for the spin averaged singlet channel. In a previous paper || 
we evaluated all of them (recently these results were confirmed, see Ref. P31 ) by taking advantage 
of conformal OPE for the process in question. The way they were obtained will be discussed in 
great detail in the next section. 

Below we outline the ideas of the second approach which allows to restore the effect of non- 
vanishing /3-function and repeat some assumptions which form the basis of this method. Its main 
idea is that the radiative corrections related to the evolution of the coupling constant represent 
the main source of large effects. In abelian theory these can be determined by resummation of 
any number of fermion vacuum polarization insertions in the gluon lines. In QCD this is no 

6 In the following consideration we will repeatedly omit the superscript Q in the LO amplitude. 
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longer true, since there is a number of other diagrams which contribute to the evolution of the 
coupling. However, this statement can be justified in the limit of an infinite number of fermion 
flavours Nf — ► oo. For moderate Nf it is clear that fermion bubbles do not produce sizable 
effects. Therefore, in a second step, which has no strong theoretical foundation, we add the sub- 
leading corrections in Nf and restore the full QCD /5-function by hand from the corresponding 



iV/-dependent coefficient |24 



A) = |iV> - \Nj - yiV c , (4) 
which is the first term in the perturbative expansion of the QCD /3-function 

A ="i =-«+««). with ¥ = s A+ (s) ,A+ -- (5) 

Thus as a first approximation one can restrict oneself to the mere summation of the vacuum 
bubble chains substituted into the NLO Feynman graphs. As we have noted at the very beginning 
in the large- Nf limit the Borel transformed perturbation series for the coefficient function T(r) = 
T n /n\(— r//3o) n can be obtained in closed form. From this the n-th order Wilson coefficient T n 
can be obtained by taking the n-th derivative (T n = l/(—/3 ) n d n /dT n \ T=Q T(T)). This technique 
gives good results for the quantities which are dominated by renormalons, although this fact can 
be traced only a posteriori by comparing the quantity in question with its exact value derived by 
some other technique. The error AT due to the asymptotic character of the perturbation series 
for Wilson coefficients can also be estimated from this Borel transformed series, i.e: 



f°° I 47T \ 

AT = ±-Im / rfrexp r T(r). (6) 

7T JO \ ttlQU / 



In the non-singlet case and in d = 4 dimensions T(r) results in a pole structure 

TM 1 jMr) ,Ri(t) R 2 (t)\ 

r(r) = fr(v + — + ^7j' (7) 

corresponding to power corrections of the form 1/Q 2 ,1/Q 4 resulting from the 1/(1 — r) and 
1/(2 — r) infrared renormalon poles, respectively. Hence R\{l)/Q 2 and i?2(2)/Q 4 yield the renor- 
malon contributions to the off forward structure functions T r . Of course, they are completely 
unphysical and should cancel against UV renormalon ambiguities of corresponding twist-4 and 
twist-6 entries due to the mixing with low-twist operators. But if one assumes that the "genuine" 
non-perturbative contribution (which is essentially related to the specific hadron) to the latter is 
small it may be used as a first approximation to the magnitude of the non-leading twist correc- 
tions. This hypothesis is referred to as "ultraviolet dominance" . In some sense this is the opposite 
assumption as the one made in the QCD sum rules approach where one takes only the phe- 
nomenological condensates into account which are accepted to be saturated by non-perturbative 
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phenomena. Power divergences due to renormalons are completely disregarded. Thus at finite 
t we can probe the non-perturbative effects seen as IR renormalons. This approach to power 
suppressed contributions cannot be claimed to be rigorous since it does not become exact in any 
limit of QCD, but rather as a sophisticated guess checked in a number of cases. In principle, 
the 1/Nf approximation kills the asymptotic freedom of QCD and thus is inadequate to describe 
the real strong interaction dynamics. However, it is by now established empirically that it allows 
to define the position of renormalons provided one replaces the QED /3-function by real /3qcd 
(it gives reliable prediction for the power r in 1/Q T corrections). Of course, this consideration 
does not introduce any new information for the processes where OPE is well established and the 
l/Q r -behaviour is known. But although OPE predictions being available for about two decades 
for DIS they have been rarely used in practice due to the absence of any systematic approaches 
for the study of non-perturbative higher- dimensional operators. In the absence of any adequate 
approach to higher twist phenomena^] it is worthwhile to take advantage of the renormalon-based 
method and try to fix the dependence of the power corrections on the momentum fraction vari- 
ables up to some numbers that fix the overall normalization. When naively applied for various 
QCD observables the absolute magnitudes R(l) were not unreasonable 0. All of this suggests 
that it makes sense to apply the renormalon estimate to our problem. We do know at least that 
for ( — > the imaginary part of our amplitude is related to the usual structure functions for which 
the renormalon analyses gave sensible results PS|. 

According to the hypothesis of UV-dominance the shape and magnitude of the higher twist 
corrections is determined by the intrinsic ambiguity of the summation of the perturbative series 
and looks like 



F r (uj, () = Jdx |T r (w, x, C, Q 2 \a s ) + 0i^A£_4(x, u, C) + 6 2 ^A 



f w _ 6 (x, W ,C) \O r (x,C), (8) 



where A 2 = fj, 2 e c (see below) and 9i are adjustable parameters which have to be fitted to data. 

Of course, the power dependence in Q 2 is defined modulo logarithms of the hard scale Q which 

are governed by the renormalization group equation. Unfortunately, for multiparticle operators 

the evolution equations are extremely complicated since the rang of the anomalous dimension 

matrix grows with the moment of the correlator involved due to an increase of the number of 

local operators mixed by renormalization. Therefore, the exact evolution equation are of the 

Faddeev-type rather then of DGLAP form. The solution of the problem at least for the non-singlet 

twist-3 sector can be found by considering the integrated quantities which depend only on one 

7 It is well known that due to the higher dimensionality of the quantities involved even the state-of-the-art QCD 
sum rules cannot claim an accuracy better than 50% — 100%, although sometimes people quote smaller errors what 
is of course not legitimate. Moreover, only few lowest moments can be studied. Lattice QCD cannot help in this 
case for the time being due to e.g. unsolved questions of operator mixing on the lattice. 
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argument, and going to the multicolour limit. In the present case the situation is complicated by 
the exclusiveness of the kinematics which results in a mixing of the operators with total derivatives 
so that even on leading twist level the anomalous dimensions turn out to be matrices not numbers 
(to say nothing about higher twists). 

In the present paper we present, apart from results on higher order corrections derived from the 
application of conformal OPE, the theoretical predictions for all-orders coefficient functions within 
the NNA approximation and a model for the momentum fraction dependence of the higher twist 
corrections coming from the assumption of the UV dominance of the non-leading twist matrix 
elements. A complete and thorough numerical analysis will be presented elsewhere together with 
the implication of different models for the non-forward distribution functions. Since the relative 
contributions of the non-perturbative and perturbative corrections we have calculated will depend 
on the functional forms assumed so will the final results. For the time being the status of such 
models is still very unsatisfactory. 

2 Beyond leading order. 

It is well known since for a long time that the conformal invariance of the theory puts severe 
restrictions on the possible form of the amplitudes. In massless field theory this is true only 
at tree level, but fails when the interaction is switched on due to renormalization effects (we 
discard for the moment the effect of gauge fixing in gauge theories since it is not of relevance 
for the physical sector). The latter can be divided in two classes: i) renormalization of the field 
operators and ii) running of the coupling constant. While the first one is not too dangerous 
since after the redefinition of the conformal representations by shifting the scale dimensions of the 
fields, given originally in terms of the canonical dimensions, by the anomalous ones, the theory 
respect conformal covariance. However, the second effect inevitably breaks conformal symmetry. 
Therefore, supposing the existence of nontrivial zero g* of the /^-function (f3(g*) = 0) conformaly 
covariant OPE can be proven to exist even for the interacting theory. Below we will shortly outline 
some of the points which are of relevance for our further discussion. 

2.1 Conformal OPE and non-forward processes. 

There exists a complete basis of twist-2 conformal operators 0, which are labeled by the conformal 
spin j and the scale dimension dj = di + d 2 +j, constructed from products of the fields of dimension 
di and spin (z/j = di + Sj — |): 

j 

CV-wy = S id H+1 . . . id m <p dl ^ C n (v u u 2 )id^ . . . id H _J d N _ n+1 ...i d N 4>d 2 - traces, (9) 
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with d =d + d and d=d — d and coefficients C n {y\,v-i). Contracting this expression with 
light-like vectors (its dual n* = p is a null vector^ along the opposite tangent to the light 
cone, defined such that p 2 = 0, np = 0) (Oji = ++ ,,, + - 3 ) we obtain immediately the well know 
expression for the conformal operators 



3l = {id+yfaP? 3,1/2 5 [d+ /d + j <fe, (10) 

where P are the usual Jacobi polynomials |26|, |28[| . 

From dimensional counting alone we can easily write for the T-product of two local (scalar) 
currents J A and J B 

°° / 1 \ \{ d A+d B -dj+j) oo 

r { Ja(o) j b (x)} = E - E • • • /*«<V^,-(o). (ii) 

Since the operators (9 transform covariantly with respect to the algebra of the conformal group 
there exists a relation between the coefficients. This relation can be found by applying, for 
instance, the generator of the special conformal transformation to both sides of Eq. dTT|) . This 
leads to a recurrence relation with the following solution: 



, . T(d j +j)T(Ud B -d A + d j -j + 2l)) 

Cn = (-1)^ 3 V2 -, — -^C.j. (12) 

r(l-j + l)T(d 3 + l)r(l(d B -d A + d 3 -j^ ' 



Performing the sum over I we get finally 



T{J A (0)J B (x)} (13) 

x ff dA O m ... w;i (0). 



v r ( M " {dA+dB ~ dj+j) r r F (\{d B -d A + d 3 + j) 



Thus the advantage of the conformal OPE for the non-forward processes is that the Wilson 
coefficients are fixed entirely by symmetry up to Cjj-coefficients which can be fixed from the 
forward matrix elements known, for instance, from DIS. 

For the case of two electromagnetic currents in order to avoid complications which arise due to 
the fact that the latter carry Lorentz indices it is enough to consider the trace and the antisymmetic 
part of the amplitude. In this way the off-diagonal part (in the conformal basis) of the amplitude is 
fixed unambiguously. (Remark: The difference between the non-diagonal analogues (Fy and Fy) 
of the forward structure functions F\ and F2 comes entirely from the forward Wilson coefficient 
function of DIS since as long as the former carry the same operator content the corrections to 



8 In what follows the plus and minus indices in the place of the Lorentz indices refer to a convolution with the 
vectors n and n*, respectively. 
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eigenf unctions are the same in both cases). Thus, we have to leading twist accuracy 



T% J M (0)J„(x) (14) 



where J 1 labels the polarization of the amplitude in question, and iF\ is a confluent hypergeometric 



function [£7], gg]. The projectors are defined as follows VY V = for the spin averaged case and 
Vt v = itfiv+- for the spin dependent one. 

In the interacting theory the conformal operators will mix under renormalization in MS-type 
schemes so that covariance will be lost. Since the mixing matrix is triangular due to Lorentz 
invariance (only the operators with the same Lorentz spin I can mix with each other) its eigen- 
values are given by the diagonal matrix elements which coincide with the anomalous dimensions 
known for forward scattering t|\ Thus, the renormalization group equation for the latter must 
be diagonalized first. This can be done by a finite transformation expressing conformal operators 
Oji of the free theory through the multiplicatively renormalizable ones denoted by Ojf. 

O sl = j2Bjkd jh (15) 

k=0 

with the transformation matrix B. It is expressed in terms of the off-diagonal 7 ND matrix elements 
of the anomalous dimension matrix 7 of the conformal operators Ojk- On the other hand 7 ND 
can be fixed entirely from the constraints coming from the algebra of dilatation and conformal 



generators, namely |29|, |3(J |) + J c (l), 7 = for (5 = (all entries will be specified below). As 
was shown in Ref. |3(J this holds true in every order of perturbation theory for (3 = 0. Therefore, 
in the interacting theory conformal OPE is of the same form as Eq. flI4|) provided we replace 
the canonical scale dimensions d c ^ n of the operators by dj = c/^ an + 2*7j" and rotate the tree level 
conformal operators to the covariant ones. 



Taking into account all that has been said above we sandwich Eq. (14) between hadronic states 



and get after some simple algebra the following result for the Fourier transformed amplitude 



i j=0 



l+j+*7j, 2 + j + l 7 j 
4 + 2j + 2*7f 



u( ){h'\ i d r jj \h), (16) 



where 2-^1 is the hypergeometric function P7| , |28| , and l Cf is the Wilson coefficient for forward 
scattering. 

Since the product of the coefficient function and the eigenfunctions of the evolution kernels is 
a scheme independent quantity we can use any scheme we want. However, for other schemes the 
manifestly conformaly covariant form of the OPE ((TBI) will be hidden. 
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2.2 NLO corrections. 

The main idea of the conformal approach for the evaluation of the coefficient functions consists in 
the combination of the information coming from conformal OPE and from ordinary factorization 
theorems. On the one hand we have the prediction for the amplitude from conformal OPE outlined 
in the previous section with an input gained from the forward scattering coefficient functions and 
anomalous dimensions. On the other hand NLO factorization tells us that the total a^-correction 
is generated by two sources: i) a s -corrections to the coefficient function and ii) a^-corrections to 
the evolution kernels which lead to a modification of the eigenfunctions of the evolution equation of 
order 0(a s ). Thus, since the latter can be fixed from conformal constraints we can get information 
about the former from the combined use of the first and second representation of the amplitude. 
Schematically, it reads 

irp i^-iCOPE irpkemel (17) 

where the conventions we have used are self-explanative. In this section we will only deal with the 
s-channel contribution. The crossed amplitude can be obtained at the end by a mere substitution, 
namely x — > ( — x. 

The general conformal decomposition of the non-forward distributions looks like 

oo 

l O(x,C)= E Y, ik <t>M^\ a s){ti\ k O j5 \h), (18) 

k=Q,Gj=0 

where the partial conformal waves are generalized, beyond one-loop level, to non-polynomial 
functions which are the subject of the constraints. 

The correction to an eigenfunction is defined completely in terms of the I?-matrix 

*=£aS) c|(2i - i)b - w=^ +i <£§w(h - < i9 > 

which in the conformal limit is defined completely in terms of the special conformal symmetry 
breaking matrix 7 C p9| , |9|] via the following relation 



B= 1 * c , with JM= Mj \ 'A fc+i, for VM- matrix, (20) 
1 + J*f a{j, k, 3/2) 

where a(j, k, v) = 2[(j + +2) — (k + l)(k + 2) + (2u - 3)(j - k)]. To obtain fa in the simplest 
way one should derive the differential equation for the latter by using the eigenvalue equation for 
the Gegenbauer polynomials as well as the identity a(j, k, 3/2)Bjk = —{ / y c B}jk which follows from 
the definition of the I?-matrix (|20|). 
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For the singlet quark case, which is relevant for our present analysis, they have the form[] 
Q %(x, CM = {^(x -y) + ^ Qk ^ y, C)} ® N ^ )(3+1 



y l 1 y 



C V C, 



Cf 2|-1 



(21) 

where 1/ = z/(fc), fc sums up the parton species (Q,G), and we introduced the shorthand notation 
® = jdy. 

To fix the normalization of the amplitude coming from conformal OPE we insert the confor- 
mal partial waves expansion into the LO result (0) which followed from factorization and get as 
prediction for the DVCS function: 



FEo(u, C) = -2 E B U + + 2)^ +1 2 F! 

3=0 



+1 p fl+j, 2+j 
4 + 2j 
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c)m Q om. (22) 



Here the conformal matrix elements of the non-forward distributions are defined in terms of the 
conformal operators sandwiched between the physical states in question 

The gluonic case differs by the replacements | — ► | and (J, I) — > (j, /) — 1, and J 1 stands now for 
the Lorentz rather than Dirac structures: 




(24) 



2.2.1 Polarized sector: gluons. 

As we have seen above the correction to the eigenfunction ® G ([> can be found with the help of the 
special conformal anomaly matrix in the quark-gluon channel. For the QG-sector the derivation 
of the conformal Ward identities can be performed in abelian gauge theory following the same 
reasoning as for the non-singlet QQ-sector. This is true since in the conformal limit, the symmetry 
breaking parts of the kernels do not contain the Casimir operator Ca of the adjoint representation 
of SU(3). Moreover, for non- vanishing (3 the correct results can be reconstructed by substituting 
the QCD /3-function for the QED one and, therefore, this treatment is sufficient to derive reliable 
NLO predictions in QCD. The special conformal anomaly matrix in LO looks like 

QG r = ^i QG iXi)}_ + Q% ] , (25) 

with b jk (l) = 6 jk {2(1 + k + 3)5 jk - [1 + (-l) j ~ k }(2k + 3)} and QG 7 being the quark-gluon mixing 



anomalous dimension matrix. Thus, the only difference compared to the analysis of Ref. [29 



The subtleties related to the spectral restriction < y < £ will be clarified in section |3^2|. 
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Figure 1: One- loop Feynman diagrams for the QG special conformal anomaly matrix. The 
empty vertex corresponds to the modified Feynman rules according to Eq. (|26|) while the full 
point denotes the usual quark-gluon vertex. The blob with cross stands for the quark conformal 
operator which mixes with the gluonic one due to renormalization. 

is that now we have to evaluate the mixing renormalization matrix Q G Z* = i^ G Zm + ... of the 
conformal operator and the operator insertion 

A g l = g— I d d x2x.C(x), (26) 



dg 

coming form the variation of the action with respect to the special conformal transformation 
§c s = _&[A£] + . . .. This is the only source of symmetry breaking, in contrast to the situation 
in the QQ-channel where subtleties arise due to the quark equation of motion. Thus, only the 
renormalization problem of the above mentioned operator has to be solved^ 

i [ Q O jl }[A s _] = i[ Q 3 A 8 -] + i E { QQ Z*} [U-i] + i E { QG Z*} [ G O kl ^\ + GVC, (27) 

fc=l 3 k=l 3 



where GVC stands for the gauge- variant counterterms |JT], which are not of relevance for the 
present discussion since they cannot affect the gauge invariant quantities we are interested in. 
To lowest order in the coupling constant we have to calculate the diagram in Fig. [l] with one 
of the vertices being replaced by the operator insertion i[A£]. To 0(a s )-accuracy we are limited 
to, it reduces to i[A£.] = igfi e f dx2x-ipflip an d thus in the Feynman rules it results in a mere 
differentiation with respect to the external gluon momenta of th e last graph in Fig. |T| but 
with the familiar quark-gluon vertices. Since the differentiation of the Gegenbauer polynomial in 
momentum space is proportional to the 6-matrix acting on the latter 

( d +4A (p + +p' + y^ ( p z ± ^4) = th k (i)(p + +p' + r i q ( p -^4) , (28) 



dp+ dp' + ) 3 \P++p' + J Z^q 3 \P++P+ 
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Square brackets mean minimally (MS) subtracted operators [[H] 
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we finally get 



Thus, the special conformal anomaly can be obtained immediately and reads: 



(29) 



(30) 



Following the discussion preceding this section we find from Eq. (|30D that the corrections to 
the eigenf unctions are completely expressed in terms of the shift operator: 



QG $(x, y, C) = (T - V)S(x, z) ® QG K A (z, y, (), 
where the generalized evolution kernel reads[^] 

Q G K A (x, x', C) = 2N f T F G{ 12 {x, x-(,x-x'), 
and S generates the shift of the Gegenbauer polynomials index: 



S{x,y)®[yyr- 1 *C»{2y-l)= ± 



[ X xY- 1 2+PC v j +p {2x- 1). 



(31) 



(32) 



(33) 



p=0 



In Eq. (|3lD X is an identity operator and T> extracts the diagonal part of any test function r(x, y) 
in its expansion with respect to a basis of Gegenbauer polynomials CJ, i.e. 



J* dxC]{2x - 1) | ^ | t(x, y) = Y. Tjk ' 1 



<%{2y - 1). 



The diagonal matrix elements of the shift operator S which will be used below are for instance 



S n = 3^(2 +j) - 2^(4 + 2j) - V(l) - 



2 + 2j 
2 + j 



(34) 



Here and below ip^ n \j) = ^^rlnr(j) denotes polygamma functions. The evolution kernel intro- 



duced in Eq. fl52"D can be diagonalized with the following identity 



c 



JdyCf -A QG K J 



(y,x,C) = QG 7f^L (2^-1] with^ 7 f 



12N f T F 



'3 3 



- (35) 



\ 7 ' 0' + i)0' + 2) 

Technically, to get this equality one has to multiply both sides by the factor (yy) 2 , differentiate 
three times with respect to y and use the relation 



{yyfC) t\{2y - 1) 



(36) 



11 In this and subsequent formulae we followed our previous convention and introduced the generalized step 



functions 



00 da_ m TT" 

-00 2-rri 



Uk=i ( ax k - i + My 
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The coefficient ® G ^ A in front of the Gegenbauer polynomial coincides with the DGLAP moments 
of the kernel G:G K A (y, x, 0) up to the common factor 6/j, which arises as a result of the conventional 

3 

definition of the Gegenbauer polynomials, namely, the coefficient of x J in CHx) is 3/j times the 

5 

coefficient of a;- 7-1 in CjL^x); an additional factor of 2 comes from the argument of the polynomial. 

Since now we know the a^-correction to the eigenfunctions of the two-loop QG-evolution kernel 
we can find coefficient function in the corresponding channel. On the one hand conformal OPE 
( |16D to NLO gives the following prediction 



G- 



NLO 



co , 

E^C? + 1 >3 + 2)^ +1 i [ G C%F 1 

2-^1 



3=0 



+ Q Gj A 



dp 



p=0 



l+j, 2+j 
4 + 2j 

1+j + p, 2 + j+p 
4 + 2j + 2p 



(37) 



^)}{h'\ G Of 3 \h), 



with the forward coefficient function 1133] 



J 



In 



-Q 2 



</>(j + l)+^(l) 



(38) 



j(j + i)(j + 2) { i? 

On the other hand factorization theorems tell us that the total a s -correction to the amplitude 
reads 

G ^nlo = g%) (w, y) ® QG ^(x, y, C) ® °O a (x, C) + G T A (u, x, C, Q 2 \a s ). (39) 
where g T,4n is the quantity in question. To proceed further we equate the right hand sides of Eqs. 



(i) 

and (|39|) and use the formulae 

CO 

T (0 ) (w, y) ® %, z) ® QG ^ A (z, a;, C) ® G O a (x, C) = -2 £ QG 7/^(j + 1, J + 2)^ +1 



(40) 



l+j, 2+j 
4 + 2j 



2^1 



p=0 



3=0 

/l+j, 2+j+p 
I 4 + 2j+2p 



and 



T (0) (w, 2/) ® ln(l - iuy) QG K A (y, x, C) ® G O a (x, () 



x\L jj2 F 1 



l+j, 2 + j 
4 + 2j 



/ dp 



2^1 



p=0 



-2j2 QG lfB{j + l,j + 2)a< 

3=0 

1+j+p, 2 + j 
4 + 2j 



(41) 



«C)K/*Wil^>- 



Here it is easy to recognize the first terms in the curly brackets in Eqs. fl40|) and (|4T|) as the 
diagonal parts of the shift operator and the logarithm Ljj = —[if)(j + 1) — if)(l)], respectively. 
Extracting only the non-diagonal part from these identities, and adding the equation which comes 
from the forward coefficient function: 

'l+j, 2 + j 
4 + 2j 



2j2 G Cf(Q 2 /p 2 )B(j + l,j + 2)uj> +1 2 F 1 



3=0 



u(\ (h'fO A jj \h) 



(42) 



T(p)(u,z) 



!® G K A (z, x, C) In -5! + [V ln(l - zoo) - 1] QG K A (z, x, C)| ® G A (x, C), 
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we can restore the form of the amplitude coming from OPE. Subtracting from the OPE result the 
correction to the eigenfunction coming from the usual factorization approach we get the coefficient 
function we are interested in 

G T(4(u;, x, C, Q 2 ) = -^T ( o)(w, v) ® [ QG K A (y, x, In ^ + [In (1 - yu) - 1] QG K A {y, x, C)}. 

11 (43) 
By subtracting the logarithm of the hard scale we obtain the net gluon coefficient function in the 
MS scheme. 



2.2.2 Polarized sector: quarks. 

In the quark sector there is no need to calculate the special conformal anomaly matrix for extended 



. It 



kinematics since it is known for the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) case 
was shown in Ref. [Q that the continuation to the whole x/C,, y/C-plane is unique. To leading 
order this results in the substitution of the ordinary ^-functions in the kernels by the generalized 



ones 



namely 



e(x-y) (x y\ o 
k{x,y) — — * k \7'7) Q n( x ~y> x ~ 0, 



(44) 



with a crossed contribution obtained by the substitutions x — > (—x, y — > (—y. The correspondence 
between extended kernels and ER-BL ones is the following K(x,y,( = 1) = — Ver-bl(^) y)- 
This means that, in principle, one could get the non-forward evolution kernels from their ER-BL 
analogues which are known for over fifteen years rather then calculate them once more. Coordinate 



space results could be reconstructed with a help of the Fourier transformation [p2 



Taking into account the effect of non-vanishing /3-function results in a modification of the 
eigenfunction of the NLO evolution kernel by an additional term ~ — j3o(I — T>)S, namely 



QQ <P(x,y,0 = (l-V){S(x,z) 



QQK A (z,y,C)-^S(z-y) 



Q G(y,x,Q , (45) 



where 



QQ K A (x,x',() 
= C F 

Q G(x,x'X) 
= C F 



X & l x (x, x -x') + - — -Q° u (x -(,x-x') + e? u (x, x - C, x - x') 
x — X x — X 



(46) 



(47) 



x'-( 



x — X 



In 



x 



e° 11 (x-(,x-x') + 



X 



In 



x — X 



x — X 



X 



Q^x.x-x') 



From conformal OPE it follows that 

oo 

-2j2 Q C A (Q 2 /^ 2 )B(j + + 2)^ +1 2 F 1 

3=0 



1+J, 2 + j 
4 + 2j 



C\(h'\ Q A \h) 



(4J 
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T(p-)(u,z) <g> \8(z - x) - 



2^ 
3 



QQ K A (z, x, C) In + V ln(l - zu;) 9K^(z, x, () 



-- Q K b (z, x,()-V QG(z, x, C) + -5(z - *)] } ® ^(x, C) 
where the forward coefficient function is taken from Ref. 

Q Cf(Q 2 /^) = 



1 + 

2tt 



25 14 (;) - 2S 2 (j) + ^(j) (~ + J— + 1 



+ 



J + l 



J + 1 J + 2 
?+f?-5iO')-5iO" + 2))lB 



(49) 



-Q 2 



Here the additional kernel reads 



Q K\x,x',Q = C F 



x-C 



X — X' 



-e° n (x -c, x -x') + -^—e° n (x, x - x') 

x — X 



(50) 



In the derivation we have taken into account that the diagonal matrix elements of the matrices 
introduced above read 



1 Q K b 

— Q G- 
C F n 

1 QQk a . . 



2[^0" + 2)-V(l)]-2, 

^ 1) (J + 1)-C(2)-[^(J + 1)-^(1)] 2 , 



ir Q V = 1>U + 3 ) + V>0' + 1) " 2V(1) " ~- 



(51) 
(52) 
(53) 



To recover the same form of the result as in Eq. ( f49|) it is useful to modify these equalities by using 
the identity 2S 1}1 (j) = Sf(j) + S 2 (j) @, where ^ m (j) = El=o S m (k) / k l , S m (k) = £*L V*™ 
ip^'ij + 1) = C(2) — S2O) HH, and by noting that with these conventions = —ICpSx^ij). 
Thus, the same steps as before give 



Q T 



CQ 2 ) = -^T {0) (u,y)®\^K A (y,x,0 In 



2tt 



jQ 2 

^ 2 



+ In (l-yu) QQ K A (y, x, C) - C)}- 



(54) 



In the limit ( = 1 we obtain from Eq. ( [54] ) the results of Refs. [[43] for the 7r — > 77 transition 



form factor F, 



7T77* 



2.2.3 Unpolarized sector: quarks and gluons. 

For the quark vector channel the only additional term comes from the forward DIS coefficient 
function |35[] and reads Cp/{j + l)(j + 2). It is easy to check that this eigenvalue corresponds to 
the evolution kernel O^n, namely 

/ dxC i 2 ( 2 f " l ) e ?u(*.* - C* - y) = {j + 1 ^ J + 2) c f [ 2 \ - l ) • ( 55 ) 
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As we have noted in the previous subsection in order to verify this and similar eigenvalue 
problems there is no need to perform the integration rather one should differentiate v + --times 
both sides of the equations multiplied by (yy) u ~^ with respect to y and use the eigenvalue equation 
for the Gegenbauer polynomials of index v. To reduce the mult i- argument G-functions to the two- 
argument ones it is enough to use the following relations 

&ijk( x u x 2 , x 3 ) = {©"-ifcO^ x 2 , x 3 ) - <3>l~\ k (x 1 ,x 2 , x 3 )\ 

= — - — {x29i- ljk (x 1 ,X2,x 3 ) - x 1 Q i y j _ lk {x u x 2 ,x z )\ . (56) 

Thus, the problem to find the off-forward analogues to some forward functions is reduced to the 
solution of the inverse problem of determining a potential from the known eigenvalues 



J dxCf^-lj K(x,y,C) = Ej 



[y)C] [2 V - - 1 j , (57) 

which in general turns out to be quite nontrivial. If this problem would be solved in general we 
would be able to reconstruct the diagonal part of the kernel in the basis of Gegenbauer polynomials 
from the corresponding forward analogue and thus the whole NLO corrections to the ER-BL-type 
evolution kernel and not only its off-diagonal part. 
Observing that 

e° m (x, x-C,x-y) = QQ K v (x, y, () - Q K b (x, y, Q - \s{x - y), (58) 

and that the corrections to the QQ-eigenfunctions in the spin-dependent and -averaged cases are 
the same we easily obtainQ: 

a sm ( . . ^ (qQt^V 7.. _ /\ i_ zS 

+ [\n(l-yco) + l] QQ K v (y,x,()- Q G(y,x,()}, (59) 
^T (0) (u, y) ® {® G K v (y, x, C) In ^ - \ r 

+ ln(l-ycu) QG K v (y,x,0}. (60) 

The calculation of the gluon coefficient function quoted here does not introduce any new specific 
features and its derivation runs along the same line as before. The momentum space evolution 
kernels Q(Q' G )K V involved above read 

QQ K v (x,x'X) = QQ K A (x,x'X), (61) 
Q G K v (x,x',() = QG K A (x, x', C) - AN f T F ( * ~ ^ Q° lu (x, x-C^x-x'). 

X [X (,J 



«T{> (lo, x, C, Q 2 ) = -gT (0) (u, y) ® [ QQ K v (y, x, () In - ^K b (y, x, () + S(y - x) 



? T ( Y } (u, x, C, Q 2 ) = -^T (0) (u, y) ® ® G K v (y, x, C) In -if- - ± \ QG K v (y, x, () + °" G K A (y, x, () 



12 



Note a misprint in Eq. (40) of Ref. || where the sign in the brackets of the fourth term should be changed. 
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2.2.4 Remark. 



To get the corrections beyond NLO one should expand the amplitude coming from conformal OPE 
( |T6D up to the required n-th order. However, these results would correspond to the special scheme 
where the conformal covariance is preserved. To be able to make the predictions in MS-type 
schemes we have to perform finite renormalization with 5-matrix which has to be known to n-th 
order. The latter requires evaluation of the special conformal anomaly 7 C in the corresponding 
channel. For instance, to NNLO this calculation goes along the same line as developed in Ref. |37j 



for the two-loop non-singlet exclusive ER-BL evolution kernel, provided we replace the ordinary 
Feynman rules, by modified ones similar to those used by us in section 2.2. 1| . (Technically this 
analysis is of the same complexity as in the paper mentioned above). However, even then, this 
prediction will be only valid for the special case of vanishing /5-function /3 — which, of course, 
may turn out to be far from reality. Therefore, it is worthwhile to study the importance of 
conformal symmetry breaking effects introduced by the running of the QCD coupling constant. 
This subject will be addressed in the next section. 



3 Resummation of fermion vacuum insertions. 

In this section we will present results for the resummation of the fermion bubble chains in the NLO 
coefficient functions of the DVCS amplitude. The idea of the NNA approximation is based on the 
observation that the effects related to the evolution of the coupling constant can be represented as a 
source of potentially large perturbative corrections. Its extraction can give important information 
on the higher order perturbative contributions. 

The resummation of the fermion loops in the NLO coefficient function leads to a factorial 
growth of the series 

oo 

T " T (o) = £4n!(-a s A)) n , (62) 

n=0 

which convergence radius is zero. Therefore, due to the asymptotic character of perturbation 
series the naive resummation is meaningless and to get reasonable predictions we should truncate 
the series at the order no at which the ratio of two successive terms is of the order 1. The first 
neglected term will give an estimate for the ambiguity of this approximation. This ambiguity is 
known to be power suppressed. This last point, discussed at length in the Introduction, will be 
used below for the construction of a model for the non-leading twist corrections. 
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(6) 



(c) (d) 



Figure 2: NLO Feynman diagrams (s-channel only) with resummed fermion vacuum insertions. 
3.1 Coefficient function. 

In the following we adhere to the convention b = — (3$ for the first coefficient of the /5-function. 
After resummation of the fermion vacuum polarization blobs (the sum is defined up to infinity in 
the sense of a principal value (PV) prescription for the poles in the Borel integral, although in 
practice the truncation of the series at its minimal term is completely equivalent to this) in the 
NLO coefficient function (see Fig. |2|) and with appropriate renormalization we get as result for the 
coefficient function for the -T-amplitude (the details of this calculation can be found in Appendix 
A) 

T r (u, x, C, Q 2 \a s ) = PV/^ *Le-*H** {T r (0, r\u, x, (, Q 2 ) - f r (-r, 0\ u , x, C, Q 2 )) , 

(63) 

where|^| T r (r, 0) is a Borel transform of T r with respect to the first argument. Note that the 
second term can be expressed in terms of the integral of the function T f (t, 0) itself 



T 



" — T r (r,0). 

T 



A straightforward calculation leads to the following representation of the Borel amplitude 

T r (e,r) 



F(e,r)LH ( e , r ), 



with a common factor F 

,*2 



F(e,r) 



-Q 2 



r(i + e )r 2 (2-e) 

T(4 - 2e) 



7TT 



T(2-e) 



sin(Trr) T(l - e + t)T(3 - e - t) 



(64) 



(65) 



(66) 



13 In what follows we omit the dependence of the Borel transform T r (e,r) on the momentum fractions as well 
as other kincmatical variables. 
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(/i* 2 = 4/T/i 2 ) and a function D which is specific for the particular channel (r = A, V). As can be 
seen form Eq. (|63|) in practice we need the function F only for vanishing first or second argument. 
We can easily get 

F(0,r) 



'/i 2 e CN 



1 



(67) 



y_ Q 2) (1 _ r)(2 _ r )- 

where C differs for different subtraction schemes: Cms = § + hi47T + ip(l), C^jg = |. 

/ x 1 (1 - t) T(4 - 2r) 

F < T - 0) = 6 (2- T ir ( i + V)P(2- T ) ' (68) 

The contributions of individual diagrams to the function D in the polarized and unpolarized 
channels are given in Appendix A. Below we just quote these expressions for those particular 
values of its arguments which enter Eq. (p3|) . 



3.1.1 Polarized case. 

D a (0,t) = T {0) (u,x] 
1 



(1 - XOj) n 



l + r 



+ (2 



l + r 

„ 1 — XUJ 



(2-r 



1 — axu 



l + r, l + r 



T 



l + r, l + r 



D a (t,0) = T {0) (u,x 

2-r 
~1 - r 



3 - r 2 2-r 



1 - r 1 - r 



l-(l-r) 



1 — XLU 



2 + r 



2 + r 



(69) 



1, 1-r 
2-r 



1 + (1 - r) 



1 — £C<J 



1, l-r 
2-r 



(z - (W 



(70) 



3.1.2 Unpolarized case. 



D y (0,r) = T (0) (u;,x 



'l-xtu) T l + r 
1 



2"(2-r(l-r)) 



1 — XOJ 



D v (r,0) = D a (t,0). 



2 + (2-r(l-r)) 



1 — XLU 



2^1 



l + r, l + r 
2 + r 



l + r, l + r 
2 + r 

(x-cv\ 

1-Cw J 



XOJ 



,(71) 
(72) 



The general feature manifested by the results we have just derived is the existence of a few IR 
renormalon poles in the Borel plane. This is a common point for all space-like processes which 
tells us what kind of power suppressed contributions should be added to make physical quantities 
free from ambiguities. 



3.2 Extended ER-BL evolution kernel. 

The question of the NNA-corrections to the eigenfunctions of the generalized Efremov-Radyushkin- 
Brodsky-Lepage evolution equation was studied by us in Ref. (see also |4T|). They are given 
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by Eq. (|19|) but now the conformal anomaly 7 C in the S-matrix should be shifted by the term 
proportional to the /^-function 7 C — > 7 C 



2-b. This follows from the conformal constraints extended 



to the case of a running coupling [29]. It was shown there that the eigenf unctions of the extended 
ER-BL evolution kernels with the fermion renormalon chains resummed to all orders leads to a 



shift of the index | of the Gegenbauer polynomial by the amount —@(g)/g: 



1 - 



x 



c 



3/2-P/g 



3 



(73) 



The NNA-anomalous dimension matrix is diagonalized in this basis 



NNA 



(74) 



and its diagonal matrix elements coincide with the forward anomalous dimension of the DGLAP 



evolution equation with fermion-loop insertions resummed to all orders jig, [39|, |4(], |4l 

(l-r) 2 (2-r) 



QQ 'NNA 
' 3 



47T 



00 



C f F(t 



,0){l 



(l-r + j)(2 
+ (4-2r)[V>(2- 



■r + j) 

r + j)-^(2-r)]}. (75) 



Since the LO evolution kernels are the same in the polarized and unpolarized cases the eigen- 
functions and eigenvalues are also the same. But, unfortunately, the eigenfunctions obtained by 
resumming only the fermion bubble chains with consequent restoration of the full /3-function has 
nothing to do with reality since as has been observed in Ref. [52"] there is significant cancelation in 
the evolution kernels between different conformal symmetry breaking parts for Nf = light flavours, 
namely, between the special conformal anomaly and the /3-function term. However, in the present 
approach while the latter is taken into account the former is discarded completely. The breaking 
of the NNA approximation for the anomalous dimension of forward deep inelastic scattering has 
been observed also in Ref. BP| which supports our conclusions. 

The above result could be seen from the expressions for the eigenfunction found in the previous 
section |2.2.2j . Adhering to the NNA we should neglect all terms in Eq. (|45D except for those ~ (3 
and to 0(a s ) accuracy (for brevity we put £ = 1) we find 



{I-{I-V)rS}®[yy]Cf{2y-l 



l- T £,3/2-T 



(22/ 



0(t 2 



(76) 



where r = f^Po, and A/}(t) = 1 + rSjj comes from the diagonal part of the shift operator which 
affects only the overall normalization. In this way we support the hypothesis of NNA since in Eq. 
( [55]) (3q stands for the full QCD /^-function rather then only for its fermionic piece up to limitation 
about its quantitative validity. 
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Note that the argument of the Gegenbauer polynomial CJ(i) is defined on the segment — 1 < 
t < 1, thus the above equalities are defined for a distribution with the support < x < ( @. 
However, if one exploits the identity 

u-l/2 

l) (77) 



1 - 



x 



c. 



X 



C 



<,l-2u 



r (I) r(i 



2z/) 



(it(tt) j+ ^ i/2 5 (i) (a-^) 



r(z/)r(j + z/ + |)r(j + i) 

and can understand the latter in the sense of a mathematical distribution, i.e. to get a meaningful 
result it should be convoluted first with some smooth function before integration over t 0, one 
can then restore the support properties to the whole range < x < 1. 



4 Renormalon model for the higher twist corrections. 



The model for higher twist contributions can be traced from the estimation of the ambiguity in 
the resummation of the perturbation series and is given by the imaginary part of the Borel integral 
which appears due to the necessity to deform the integration path in the complex plane of the 
r-parameter to escape from singularities (this leads to the undetermined overall sign owing to 
different possibilities to close the integration contour around the pole) 

2C F 1 , r°°dr 



l - Im r ^e-^ b >T r (0, r| W , x, (, Q 2 ) = ± £ 

71 JO T JTT, 



y e CN(n-2)/2 



tw— n 



-Q 2 



A^ n (x,u,0- (78) 



There are only two poles in the Borel plane which are located in the points r = 1,2. Thus, the 
ambiguities generate only two types of power corrections corresponding to twist-4 and 6 operators. 
Of course, one can generate even higher power correction by considering the perturbative diagrams 
beyond one-loop order. However, in view of the success of the leading IR renormalon model we 
will limit ourselves to the one- loop approximation. From Eqs. (|69D, ( fTTD of the preceding section 
we get for the spin- dependent sector 

AC F T( )(u;, x) ( (1 — C^)(l — xoj + (u) 



A? w _ a (x,uj,() 



1 

XUJ 



At 



{x,u,0 



b ( 

1 — XUJ 

2C F 



x - Quj I (uj(x - Quj 
(2x - ()uj(1 - (ujY 



ln(l - (to) 



(79) 



xuj(x — Quj 
1 - 



b T ^ X \(x-C)^ 



1 — XUJ 



2(1 



XUJ 



ln(l — xuj) + 

ln(l-Cw)- 
xuj{1 — (uj) 



1 



(XUJ 



XUJ 



+ 



[x — C) 3 ^ 3 



1-(UJ 



XUJ 



ln(l 



XUJ 



(XUJ 



XUJ 



(x - C) 2 ^ 2 
(x - C) 2 uj 2 \ { 



(80) 



1 - Cu)(2x - C)u - 



XUJ 
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and for the spin-averaged one 



A 



Af w _ 4 (x,^,C), 
2C F 



T, 



1 



1 + 2 



XUJ 



1 — XUj\ 1 — Quj 
(uj , 



(x — () 3 UJ 3 

1 — XUJ 



[XUJ 



(x - C) 3 o; 3 



1 + 2- 



ln(l - (uj) 
ln(l 



*1) 



XUJj 



2{i-xuy 



[XUJ 



1 — XLU I 

XUJ \ 

1-Cuj 



2 - 3xu;)\ 
xuj^uj y 
1 — xa/ 



(82) 



[x 



O 2 ^ 



Cuj 



1 — 3 xuj + 



- CV, 



We should note that the twist-4 result can be considered as more reliable since the only contribution 
to it comes from the IR renormalons in the coefficient function for the twist-2 operators. For the 
twist-6 part there is additional input from IR ambiguities in the point r = 1 from the coefficient 
function of the twist-4 correlators. Thus the above presented results for tw-6 include only part 
of the total contribution. However, since the latter is suppressed by extra two powers of the 
momentum transfer its particular shape is not of relevance for us and we have cited them just for 
completeness. 

There is another source of power suppressed contributions which comes from the non-conver- 
gent behaviour of the perturbation theory and is due to instantons |45l . However, they lead to 
strongly suppressed contributions due to high values of r and thus can safely be discarded. 



5 Summary and outlook. 

In this paper we have presented the theoretical framework for studying higher order and higher 
twist corrections based on the combination of two different formalisms. The first one is based 
on restrictions for the amplitudes of the massless theory coming from the algebra of the collinear 
conformal group in the hypothetical limit of vanishing /3-function. We restore the dependence on 
the latter with the second approach which resums the fermion vacuum insertions to all orders in 
the coupling. Taken alone this approximation is insufficient since fermion loops do not dominate 
the radiative corrections to the amplitudes. However, supplied with the idea of NNA one may hope 
that they do and that thus all important perturbative corrections are taken into account. The 
validity of this approximation can be checked only by comparison with exact results. It turns out 
that not all quantities are approximated with good quality by this methods, but only those which 
are dominated by the renormalon poles. Due to the latter the resummed amplitude is plagued 
by the uncertainties that manifest the asymptotic character of the perturbative series. It was 
established that they are power suppressed and thus mimic power corrections to the amplitude. 
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This tells us that both higher orders and higher twist should be treated simultaneously to escape 
from ambiguities intrinsic to the perturbative series. By accepting that the non-perturbative higher 
dimensional operators are dominated by their UV renormalon poles it is possible to construct rough 
model for the momentum fraction dependence of the power suppressed contributions. 

In conclusion we have resummed fermion vacuum polarization bubbles in the coefficient func- 
tion of the DVCS amplitude. It is known that NNA approximation overestimates radiative correc- 
tions. A more realistic estimate is given by the semi-sum of the exact NLO (section H) and NNA 
(section |3.1|) results. By using the UV dominance hypothesis we give an estimate of the shape 
and magnitude of the higher twist contributions. The numerical analysis will be performed in a 
separate publication. For that we have to accept one of the models for leading twist non-forward 



distribution. Several of them are already on a market [46, 47 



Note added: After the present study was finished we learned about Refs. [49], [5(| where the 
leading twist factorization for DVCS was proved to all orders of perturbation theory. 

We wish to thank M. Maul for the collaboration at an early stage of the work and A.V. 
Radyushkin for comment. A.B. was supported by the Alexander von Humboldt Foundation and 
partially by Russian Foundation for Fundamental Research, grant N 96-02-17631. 



A Fermion bubbles in the DVCS coefficient function. 
A.l Spin-dependent scattering. 

In the treatment of the polarized sector we resolved the 75-ambiguity following Braaten's recipe 
[ [|?|] , namely, we used the anticommutativity property of 75 only in the box- type diagram and 
contract the string of 7-matrices through 75, while in all other cases we did it in the other direction 
such that no commutation with the chiral matrix occurred. 

We have calculated the one-loop Feynman graphs represented in Fig. |2| in a d = 4 — 2e 
dimensional space with the following effective gluon propagator in the Landau gauge (/i* 2 = Aufi 2 ): 

<-'«>™(5)>-^)- ^ 

A simple calculation leads to the results (since the sum of all diagrams forms a gauge invariant 
set, only the g^-part of propagator has been used in the real calculation) 
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While the calculation of the self-energy and vertex-type corrections is straightforward we should 
mention some technical details about the box-type graph. As quoted above a straightforward 
calculation of the latter leads to a result in terms of the Appel function F\ p?| , However, 
the latter can be reduced to the difference of two hypergeo metric functions 2F1. The derivation 
proceeds in a very simple manner. One evaluates the momentum integral in <i-dimensions with the 
integrand given by the product of the denominators of the box diagram and the factor (pk). There 
are two possibilities to compute this integral: on the one hand we can join all four denominators 
via Feynman parameters and get at the end the result in terms of Fi. The second way is to 
represent the factor (kp) as a difference of two propagators, namely 



(kp) = i {[{x - ()p - k] 2 - [xp - k] 2 } . 



(A3) 



Then, the final result looks like a difference of the vertex-type graphs. Thus, we end up with the 
following relation between the Appel function and hyper geometric ones which we have failed to 



find in any textbook on special functions fl27| , |28 
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We put a = ne, with n being the number of the fermion bubbles inserted in the gluon line, and 
multiply the above expressions with a factor corresponding to the product of the fermion vacuum 
polarization blobs 
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Performing trivial subtraction of sub- and overall divergences and resummation, which is partic- 
ularly easy as only the end terms in the sum survive [44|, along the line of Ref. J7| we obtain the 
result given in the main text in Eq. (^) with the functions Z), corresponding to the contributions 
of particular diagrams 

1 -e 
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Although our expression for the box-type diagram differs from the one calculated (for £ = 1) in 
Ref. [41] by the factor ~ er obviously this contribution has no impact, neither on the final answer 



fl63|) since it vanishes in both limits e = r = 0, nor on the one-loop results derived in section |2.2.2 
as it proportional to e 2 . 



A. 2 Spin- averaged scattering. 

For the spin-averaged amplitude there is no difficulty due to 75 and the calculation is straightfor- 
ward 
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where the function /(er, e) is defined by Eq. ( A. 2 ). Some simple manipulations give finally 



Dt(e,r) = D v a (e,r), D*(e, r) = (e, r), D c A (e, r) = D v c (e, r), 
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